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An  Evaluation  of  an  Advection  Fog 
Prediction  Model 

I.  IYTROIHCTION 

In  response  to  Geophysical  Requirement  GR  3-76,  the  Air  Force  Geophysics 
Laboratory  (AFGL)  has  been  involved  in  research  to  develop  short-range  advection 
fog  prediction  techniques.  This  research  has  been  a  two-pronged  study.  The  first 
part  has  concentrated  on  collecting  various  types  of  data  at  and  near  the  Weather 
Test  Facility  at  Otis  ANGB,  Mass.,  to  gain  a  better  understanding  of  the  micro¬ 
physical  and  dynamic  processes  in  fog.  The  results  of  this  work  have  been  re- 
1  2  3  4  5 

ported  by  Kunkel,  ’  *  '  Mack,  and  Hanley  et  al.  The  second  part  of  this  pro¬ 

ject  was  to  select  an  advection  fog  prediction  model  and  to  test  its  accuracy  by  com- 

(Received  for  publication  5  June  1984) 

1.  Kunkel,  B.A.  (1981)  Comparison  of  Fog  Drop  Size  Spectra  Measured  by  Light 

Scattering  and  Impaction  Techniques,  AFGL-TR-81-0049,  AD  A100252. 

2.  Kunkel,  B.A.  (1982)  Microphysical  Properties  of  Fog  at  Otis  AFB,  AFGL-TR- 

82-0026,  ADA119SZS: 

3.  Kunkel,  B.A.  (1984)  Parameterization  of  droplet  terminal  velocity  and  extinc¬ 

tion  coefficient  in  fog  models,  J.  Climate  Appl,  Meteorol.  33:80-89. 

4.  Mack,  E.J.,  Wattle,  B.J.,  Rogers,  C.W.,  and  Pilie,  R.J.  (1980)  Fog  Char¬ 

acteristics  at  Otis  AFB.  MA.,  Calspan  Corporation  Final  Report,  AFGL- 
TR-80-0340,  AD  A095358. 

5.  Hanley,  J.T.,  Mack,  E.J.,  Wattle,  B.J.,  and  Kile,  J.N.  (1983)  The  Feasibi¬ 

lity  of  Using  Precursor  Aerosol  Parameters  to  Predict  Visibility  in  Fog  and 
Haze,  Calspan  Corporation  Final  Report,  AFGL-TR-83-0195,  ADA137141. 


paring  the  model's  results  with  observed  case  studies.  Then,  based  on  the  results 
of  this  comparison,  methods  to  improve  the  model  would  also  be  sought. 

Historically,  researchers  first  had  the  facilities  to  develop  realistic  numeri¬ 
cal  models  of  the  atmosphere  in  the  late  1950s  and  early  1960s.  These  initial  ef- 

C  17  Q 

forts  focused  on  many  different  areas.  Blair  et  al,  Lilly,  Ogura,  and  others 

q 

were  simulating  buoyant  convection.  Ogura  and  Charney  developed  a  model  of  a 

squall  line,  and  Orville*0  was  investigating  mountain  upslope  winds.  In  the  area  of 

1112 

the  atmospheric  boundary  layer,  Estoque  ‘  simulated  diurnal  temperature  va- 

13 

riations,  while  Fisher  and  Caplan  predicted  the  development  and  dissipation  of 

14 

fog.  A  classic  paper  by  Lilly  discussed  models  of  cloud-topped  mixed  layers. 
Since  these  primitive  beginnings,  many  researchers  using  hundreds  of  models  have 
studied  the  atmospheric  boundary  layer  and  fog/stratus  systems.  These  research¬ 
ers  have  tested  and  examined  such  items  as  higher  order  turbulence  closure  tech¬ 
niques,  methods  to  handle  the  entrainment  of  dry  air  from  above  an  inversion, 
short-  and  long-wave  radiation  effects,  three-dimensional  models,  moisture  in¬ 
puts,  subsidence,  exchange  coefficients  based  on  model  variables,  and  surface 
boundary  layer  fluxes.  However,  with  the  increased  complexity  of  certain  models 
and  with  the  addition  of  a  third  dimension,  much  more  computer  memory  and  time 
was  necessary.  With  the  technical  advances  in  computer  hardware,  however,  the 
cost  of  this  larger  computing  power  has  decreased. 

In  this  research,  the  selection  of  the  model  to  be  used  was  based  on  three  re¬ 
quirements:  (1)  the  memory  requirement,  which  could  not  exceed  that  of  a  mini¬ 
computer;  (2)  atmospheric  dynamics,  which  must  be  sufficiently  sophisticated  to 

6.  Blair,  A.,  Metropolis,  N. ,  von  Neumann,  J.,  Taub,  A.H. ,  and  Tsingou,  M. 

(1959)  The  evolution  of  a  convective  element:  A  numerical  calculation,  in 
The  Atmosphere  and  the  Sea  in  Motion.  Rockefeller  Institute  Press,  New 
York,  pp.  425-439. 

7.  Lilly,  D.K.  (1962)  On  the  numerical  simulation  of  buoyant  convection.  Tellus. 

14:148-172. 

8.  Ogura,  Y.  (1962)  Convection  of  isolated  masses  of  a  buoyant  fluid:  A  numeri¬ 

cal  calculation.  J.  Atmos.  Sci.  19:492-502. 

9.  Ogura,  Y.,  and  Charney,  J.F,  (1960)  A  numerical  model  of  thermal  convec¬ 

tion  in  the  atmosphere,  Proc.  Intern,  Symp.  Numerical  Weather  Prediction. 
Tokyo,  Japan  Meteorol.  Soc.,  431-451. 

10.  Orville,  H. D.  (1964)  On  mountain  upslope  winds,  J.  Atmos.  Sci,  21:622-633. 

11.  Estoque,  M.A.  (1959)  A  preliminary  report  on  a  boundary  layer  numerical  ex¬ 

periment,  GRD  Research  Notes,  AFCRC  (ARDC)  20:1-29. 

***** 

12.  Estoque,  M.A,  (1963)  A  numerical  model  of  the  atmospheric  boundary  layer, 

J,  Geophys,  Res.  68:1103-1113. 

13.  Fisher,  E.L.,  and  Caplan,  P.  (1963)  An  experiment  in  numerical  predictions 

of  fog  and  stratus,  J,  Atmos.  Sci.  20:425-431. 

14.  Lilly,  D.  (1968)  Models  of  cloud-topped  mixed  layers  under  a  strong  inver¬ 

sion,  Quart.  J.  Roy.  Meteorol.  Soc.  94:292-309. 


produce  realistic  results;  and  (3)  availability  of  the  computer  code.  The  first  re¬ 
quirement  was  necessary  because  the  model,  after  development,  was  intended  to 
be  tested  in  AFGL's  Man-computer  Interactive  Display  Access  System  (McIDAS), 
a  test  bed  for  Air  Weather  Service  planned  Automated  Weather  Distribution  System 
(AWDS). 

15 

The  model  selected  was  one  developed  by  Barker  at  the  Naval  Environmental 
Prediction  Research  Facility  (NEPRF).  Section  2  of  this  paper  briefly  describes 
the  model  and  a  few  programming  errors  discovered  in  the  review  of  the  model. 
Section  ?  details  the  sensitivity  studies  done  on  Barker's  basic  model  and  the  com¬ 
parison  with  his  results.  Section  4  describes  the  results  of  a  comparison  between 
the  model  and  six  case  studies  that  Calspan  Advanced  Technology  Center  developed 
for  the  Naval  Air  System  Command.  The  suggestions  and  conclusions  for  this 
model  are  listed  in  Section  5. 

2.  MODKI.  DESCRIPTION 

The  model  consists  of  three  distinct  levels  as  shown  in  Figure  1:  the  surface 
boundary  layer,  the  Ekman  Layer,  and  the  free  atmosphere.  The  surface  bound- 


Figure  1.  Layers  Depicted  in  Model 


15.  Barker,  E.H.  (1977)  A  maritime  boundary-layer  model  for  the  prediction  of 
fog.  Boundary -Layer  Meteorol.  11:267-294.. 


ary  layer  (SBL)  extends  from  the  surface  roughness  height,  zq,  to  a  height  of  2  m. 

In  this  layer,  the  fluxes  of  heat,  momentum,  and  moisture  are  assumed  constant. 

In  addition,  these  fluxes  are  obtained  from  the  flux  gradient  similarity  relations 

1 6 

that  permit  the  fluxes  to  be  described  from  known  mean  profiles  (Businger  ). 

The  Ekman  layer,  or  planetary  boundary  layer  (PBL),  extends  from  the  top  of  the 
SBL  (2  m)  to  the  base  of  the  free  atmosphere.  In  this  model,  the  top  of  the  Ekman 
layer  for  unstable  conditions  is  defined  as  the  base  of  the  inversion  capping  the 
PBL.  For  a  stable  layer,  the  top  is  empirically  determined.  It  is  assumed  that 
the  PBL  is  in  a  nonsteady  state,  is  over  homogeneous  terrain,  remains  incom¬ 
pressible,  and  has  nondivergent  flow.  At  its  top,  the  turbulent  friction  gradually 
diminishes  to  zero.  The  PBL  layer  is  topped  by  the  free  atmosphere,  which  is  in 
geostrophic  balance. 

Only  the  basic  equations  for  the  mathematical  relations  for  the  model  are  pro- 

15  17 

vided  here.  For  a  more  detailed  description,  see  Barker  '  and  Barker  and 
Baxter. 

2.1  Model  Equation* 

The  basic  equations  to  be  used  in  this  model  are  Eq.  (1),  the  horizontal  equa¬ 
tion  of  motion;  Eq.  (2),  equation  of  continuity;  Eq.  (3),  a  thermodynamic  equation; 
and  Eqs.  (4)  and  (5),  conservation  of  water  vapor  and  liquid  water  equations. 


if  ■ 1 '  -  ’>  ♦  h  <»  V 


|u  +  |w  =  0 

dx  dz 


3  (y  _  ^  o  ^  a  i  d  /  »  *  \  .  R 
Ft  =-ua7  -w3F+^(<7'w,)  +  p 


13  =  -u  Is  -  w  ^-+  £-  (FFo  -  — 

at  dx  az  az  4  p 


16.  Businger,  J.A.  (1973)  Turbulent  transfer  in  the  atmospheric  surface  layer. 

Workshop  on  Micrometeorology,  D.A.  Haugan,  Ed.,  Am.  Meteorol.  Soc. 

17.  Barker,  E.H.  (1975)  A  Maritime  Boundary  Layer  Model  for  the  Prediction  of 

Fog,  ENVPREDRSCHFAC  Technical  Paper  No.  4-75. 

18.  Barker,  E.H. ,  and  Baxter,  T.L.  (1975)  A  note  on  the  computations  of  atmos¬ 

pheric  surface  layer  fluxes  for  use  in  numerical  modeling,  J.  Appl. 
Meteorol.  14:620-622. 
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(5) 


In  Eqs.  (1)  to  (6),  the  definitions  of  the  symbols  are: 

Cp  specific  heat  during  isobaric  processes 
C  rate  of  condensation 
f  coriolis  parameter 

g  acceleration  due  to  gravity 

i  square  root  of  -1 

l  mixing  ratio  of  liquid  water 

L  coefficient  of  latent  heating 
q  specific  humidity 
R  radiational  heating  rate 
t  time 

T  temperature 

u  x-component  of  velocity 

V  velocity  vector  on  complex  plane 

V  geostropic  velocity  vector  on  complex  plane 
Vt  terminal  velocity 

w  vertical  velocity 
x  horizontal  coordinate 
z  vertical  coordinate 
p  air  density 

The  formulation  for  the  surface  boundary  layer  fluxes  is  fashioned  after 

19  20 

Businger  et  al,  and  eddy  exchange  coefficients  derived  by  O'Brien  were  used 

for  turbulent  mixing  in  the  PBL.  The  depth  of  the  PBL  during  unstable  conditions 


19.  Businger,  J.A.,  Wyngaard,  J.C.,  Isumi,  Y. ,  and  Bradley,  E.F.  (1971) 

Flux-profile  relationships  in  the  atmospheric  surface  layer,  J.  Atmos.  Sci. 
28:181-189. 

20.  O'Brien,  J.  J.  (1970)  A  note  on  the  vertical  structure  of  the  eddy  exchange 

coefficient  in  the  planetary  boundary  layer,  J.  Atmos.  Sci.  27:1213-1315. 
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is  an  adaptation  of  Randall  and  Arakawa.  For  stable  conditions,  the  depth  is  de- 

22 

termined  empirically  after  Clarke. 

The  radiation  model  is  a  modified  version  of  the  radiation  model  developed  by 
23 

Atwater  which  describes  the  upward  and  downward  diffuse  radiation  flux  at  some 
reference  level.  In  this  formulation,  the  transmissivity  effects  of  clouds,  water 
vapor,  and  carbon  dioxide  is  calculated,  and  a  cooling /heating  rate  is  found  from 
the  net  flux  divergence. 

2.2  Numerical  Considerations 

The  numerical  differencing  schemes  used  in  this  model  are  the  Crank-Nichol- 
son  scheme  for  vertical  transport  and  the  momentum  equation,  upstream  differ¬ 
encing  for  horizontal  advection,  one-sided  differencing  at  the  boundaries,  and  ex¬ 
plicit  centered  differencing  for  radiational  cooling. 

The  grid  spacing  in  the  horizontal  is  10  km;  in  the  vertical,  the  spacing  ex¬ 
pands  with  height.  There  are  20  vertical  columns  and  28  levels.  The  time  step 
used  to  compute  the  vertical  transport  has  a  maximum  value  of  2  minutes,  but  al¬ 
lows  the  time  step  for  radiational  cooling  and  horizontal  transport  to  be  10  times 
larger.  This  was  designed  to  take  advantage  of  the  largest  time  step  possible  in 
each  process. 

The  initial  values  for  V  are  assumed  steady  state  and  neutrally  stable,  and  are 

24 

calculated  as  described  by  Kraus.  This  prevents  oscillating  solutions. 

2.3  Coding  Error!* 

In  the  analysis  of  AFGL's  copy  of  this  model  and  its  coding,  a  number  of  er¬ 
rors  were  discovered.  It  was  confirmed  with  Barker  (private  correspondence)  that 
these  errors  were  also  in  the  original  model.  The  effects  of  these  errors  will  be 
discussed  in  Section  3. 


21.  Randall,  D.,  and  Arakawa,  A.  (1974)  A  Parameterization  of  the  Planetary 

Boundary  Layer  for  Numerical  Models  of  the  Atmosphere  (unpublished). 

22.  Clarke,  R.  H.  (1970)  Recommended  methods  for  the  treatment  of  the  boundary 

layer  in  numerical  models,  Australian  Meteorol.  Magazine.  18:51-71. 

23.  Atwater,  M.A.  (1970)  Investigation  of  the  Radiation  Balance  for  Blluted  Lay¬ 

ers  of  the  Urban  Environment,  Ph.D.  thesis.  New  York  University,  New 
York. 

24.  Kraus,  E.B.  (1972)  Atmosphere -Ocean  Interaction,  Clarendon  Press,  Oxford. 
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3.  MODEL  COMPARISON  AND  SENSITIVITY 

After  a  thorough  investigation  of  the  model,  we  attempted  to  duplicate  Barker's 
15 

results  to  ensure  that  we  had  a  current  version  of  the  model  and  that  it  could  re¬ 
produce  his  findings. 

Figure  2  is  the  temperature  profile  at  three  locations  that  Barker  reported  for 


Figure  2.  Temperature  Profile  of  "Typical  Case"  at  Three 
Locations  (After  Barker^) 


15 


/V, 


W*.‘!  vIVL  v 


his  "typical  case. "  Figure  3  depicts  the  liquid  water  content  (LWC)  profiles  in  the 
resulting  fog.  In  this  typical  case,  warm,  moist  air  was  advected  over  cool  water. 
Barker  specified  a  -0.04°Ckm  1  sea  surface  temperature  gradient,  a  geostrophic 
wind  speed  of  4  m  s  1  in  the  x-direction,  and  temperature  and  humidity  profiles  at 
the  upstream  boundary  that  were  held  constant.  Figures  2  and  3  show  that  the  at¬ 
mosphere  cools  below  200  m  and  fog  forms  near  the  sea  surface  and  slowly  deepens 
downstream. 

To  duplicate  these  results,  AFGL's  version  of  Barker's  model  was  run,  first 
in  its  original  form  (with  coding  errors)  and  then  with  the  corrected  coding  errors 
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Figure  3.  Liquid  Water  Content  (LWC)  Values  Resulting  From  "Typical  Case" 
(After  Barker’5) 


very  closely  with  Barker's  version  (Figure  2)  except  at  x  =  190  km,  where  the  tem¬ 
peratures  below  200  m  are  slightly  warmer  (less  than  1°)  in  Barker's  version. 

With  the  coding  errors  corrected,  we  note  very  little  change  at  x  ■  90  km  from 
AFGL's  version  with  the  coding  errors.  However,  at  x  *  190km,  we  see  slightly 
warmer  temperatures  below  25  m  and  significantly  colder  (2°)  temperatures  above 
that  level  compared  with  the  AFGL  version  with  the  coding  errors. 

For  the  formation  of  fog,  a  comparison  of  Figures  3  and  5  points  out  the  dif¬ 
ferences.  The  AFGL  model  in  the  original  form  developed  fog  at  the  same  location 
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Figure  5.  Liquid  Water  Content  (LWC)  Values  Resulting  From  "Typical  Case, 
AFGL  Version 
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(110  km)  as  Barker  reported;  however,  the  vertical  extent  was  lower  than  in  Bar¬ 
ker’s  version.  The  opposite  is  true  when  we  compare  LWCs  greater  than  25 g  kg  *. 
The  starting  point  for  both  is  the  same,  but  the  vertical  extent  of  the  larger  LWCs 
is  slightly  greater  for  the  AFGL  original  form  than  for  Barker's.  The  AFGL  cor¬ 
rected  model  output  shows  the  fog  developing  slightly  farther  downstream  and  a 
much  larger  vertical  extent  of  the  larger  values  of  LWC  than  in  either  of  the  origi¬ 
nal  versions. 

Based  on  the  above  comparison  and  the  small  differences  between  the  AFGL 
original  form  model  and  the  data  Barker  reported,  it  was  assumed  that  AFGL's 
original  version  of  the  model  was  very  similar  to  Barker's.  f 

After  an  analysis  of  the  coding  errors  that  were  discovered,  we  determined 
that  the  error  of  primary  importance  was  in  the  calculation  of  the  delta  virtual  po¬ 
tential  temperature,  A  6^.  The  model  calculates  the  bulk  Richardson  Number  from 
the  A0V‘  Then  the  model  uses  the  bulk  Richardson  Number  to  determine  the  sur¬ 
face  boundary  layer  fluxes  and  the  eddy  exchange  coefficients.  When  we  corrected 
this  error,  the  surface  boundary  layer  fluxes  increased  slightly,  and  the  eddy  ex¬ 
change  coefficients  increased  moderately.  The  result  was  to  increase  the  diffusion 
of  cold,  moist  air  near  the  water  surface  to  greater  heights.  Therefore,  as  shown 
in  Figures  4  and  5,  the  corrected  model  had  colder  temperatures  above  25  m, 
greater  vertical  extent  to  the  fog,  and  higher  LWCs  than  either  of  the  original  ver¬ 
sions. 

In  another  experiment.  Barker  considered  the  atmosphere  to  be  adiabatic 
above  200  m  while  all  other  conditions  remained  as  they  were  in  the  typical  case. 
This  greatly  increased  the  amount  of  net  radiation  leaving  the  boundary  layer.  Fig¬ 
ures  6  and  7  (after  Barker)  display  the  results  of  this  experiment.  A  comparison 
of  this  case  with  the  typical  case  (Figures  2  and  3)  shows  that  the  fog  is  more  ex¬ 
tensive,  the  LWCs  are  higher,  and  the  temperatures  below  200  m  at  x  =  190  km  are 
considerably  colder  than  in  the  typical  case. 

The  output  from  the  AFGL  original  and  corrected  models  for  this  experiment 
is  shown  in  Figures  8  and  9.  Again,  we  see  a  close  correlation  between  Barker's 
data  and  the  original  form  data  (see  Figures  6  and  7).  However,  this  time  we  see 
very  little  difference  between  the  AFGL  original  and  corrected  form  data.  Because 
of  the  expanding  vertical  grid,  the  errors  appear  large  above  80  m,  but,  in  fact, 
represent  a  difference  of  only  one  grid  interval.  In  this  case,  the  cooling  was  very 
large  because  of  the  net  radiation  loss  near  the  fog's  top.  It  completely  dominated 
the  diffusion  effects.  Therefore,  the  errors  introduced  by  the  coding  problem  were 

minimal  in  comparison  with  the  overall  results. 

17 

Barker  described  a  number  of  other  experiments  he  performed.  These  in¬ 
volved  studying  the  model's  sensitivity  to  moisture,  wind  speed,  and  sea  surface 
temperature  gradient.  Each  experiment  he  listed  was  tested  both  in  the  AFGL  ori- 
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Figure  7.  Liquid  Water  Content  (LWC)  Values  Resulting  From  Increased 
Radiant  Heat  Loss  Case  (After  Barkerl5) 


ginal  and  the  corrected  models.  The  original  model  results  compared  quite  well 
with  Barker's  results  for  each  of  these.  For  the  corrected  version,  results  simi¬ 
lar  to  those  discussed  in  this  section  were  found.  When  turbulent  diffusion  effects 
played  an  important  part,  large  differences  occurred  between  the  original  and  cor¬ 
rected  forms.  When  diffusion  was  dominated  by  other  effects,  then  both  versions 
closely  agreed. 


1.  RESILTS 


When  AFGL  initiated  its  research  into  short-range  fog  prediction  techniques 


Figure  9.  Liquid  Water  Content  (LWC)  Values  Resulting  From  Increased  Radiant 
Heat  Loss  Case,  AFGL  Version 


the  plan  was  to  test  the  effectiveness  of  the  selected  two-dimensional  fog  prediction 
model  with  data  obtained  at  Otis  ANGB,  Mass,  AFGL1'  3  and  Calspan^'  5  under 
contract  to  AFGL  collected  the  necessary  surface  data  at  and  upwind  of  Otis  ANGB. 
However,  fine -scale  vertical  sounding  data  did  not  become  available  until  late  in 
the  measurement  program,  when  there  were  no  suitable  fog  cases. 

At  approximately  the  same  time,  under  a  2-year  contract  (No.  N00019-81-C- 
0102)  with  the  Naval  Air  Systems  Command,  Calspan  Advanced  Technology  Center 
was  assembling  information  derived  from  recent  marine  fog  observations  and  anal- 


yzing  various  models  for  marine  fog  and  the  marine  boundary  layer.  Their  goal 
was  to  develop  operationally  useful  forecast  concepts.  Under  Task  1  of  this  effort, 
Calspan  provided  six  data  sets  with  which  various  numerical  model  investigators 
could  exercise  their  models.  Table  1  is  a  summary  of  those  six  data  sets. 

These  data  sets  provided  the  necessary  vertical  and  horizontal  meteorological 
variables  to  conduct  realistic  tests  on  the  AFGL's  selected  two-dimensional  model. 
In  addition,  the  six  cases  covered  a  variety  of  situations  that  enabled  strong  and 
weak  points  of  the  model  to  be  identified.  The  model's  temperature  and  dew  point 
output  results  from  each  of  these  cases  were  plotted  against  the  initial  and  verifi¬ 
cation  profiles  and  are  presented  in  Appendix  A.  A  summary  of  the  initial  data  and 
the  model  output  and  verification  data  is  presented  in  Table  2. 

The  outcome  of  the  comparison  between  the  verification  data  and  the  model 
data  was  very  disappointing.  The  model  did  not  develop  one  episode  of  fog  or  low¬ 
ering  stratus  in  any  of  the  six  cases.  In  Cases  1,  3,  4,  and  6,  where  a  cloud  layer 
was  initially  specified,  the  model  rapidly  evaporated  the  liquid  water,  produced 
100  percent  relative  humidities  below  the  temperature  inversions,  and  then  changed 
very  little  for  the  remaining  time  period.  The  model  was  then  examined  to  seek 
reasons  for  the  poor  performance.  Four  areas  were  found  that  adversely  affected 
the  results. 

1. 1  >cl  Flux  Pnergenee  of  Infrared  Radiation 

Case  6  is  an  evening  situation  where  initially  stratus  is  reported  from  100- 
200  m.  After  3  hours,  the  entire  layer  below  200  m  cools  and  the  depth  of  the  fog 
is  observed  from  the  surface  to  200m.  In  Figures  All  and  A12,  we  can  see  that, 
instead  of  cooling,  the  model  warmed  the  entire  layer  and  evaporated  the  existing 
liquid  water  cloud  droplets. 

The  procedure  used  in  the  model  for  handling  the  net  flux  divergence  of  the  in¬ 
frared  radiation  was  examined  to  find  an  explanation  for  this  occurrence.  This 
flux  divergence  is  directly  proportional  to  the  cooling /heating  rate  of  the  atmos¬ 
phere.  For  Case  6,  the  cooling  rate  for  the  cloud  top  vertical  grid  interval  (36  m 
thick)  was  calculated.  Instead  of  finding  a  positive  value  (cooling),  the  first  20-min 
period  produced  a  cooling  rate  of  -1.75°C/h  (net  warming).  In  the  next  lower  ver¬ 
tical  grid  interval  (30  m  thick),  which  was  completely  filled  by  cloud,  the  cooling 
rate  was  -0.7°C/h. 

To  verify  this  result,  temperature  and  flux  profiles  from  Flight  8  on  12  May 
25 

1976  (Stephens  et  al  )  were  used  as  input  into  the  model.  In  the  cloud  top  vertical 
grid  interval  (cloud  80  m  thick  in  interval),  the  rate  was  -0.  5°C/h.  In  the  grid  in¬ 
terval  below  this  (100  m  thick),  which  was  completely  filled  by  cloud,  the  cooling 
rate  was  -1.  l°C/h.  Although  Stephens  et  al  did  not  report  the  specific  cooling  rate 


for  this  flight,  on  similar  flights  the  upper  two  layers  had  IR  cooling  rates  of  1.  0 
to  2.0°C/ h  respectively. 

The  Calspan  results  and  the  results  from  the  Stephens  et  al  data  showed  that 
the  magnitude  of  the  cooling  was  correct  in  the  model  but  the  sign  was  incorrect. 
The  method  used  in  the  model  for  handling  the  IR  radiation  process  was  examined 
in  detail,  but  no  obvious  procedural  errors  could  be  found  to  xplain  the  results. 

At  one  point,  it  appeared  that  there  was  an  incorrect  sign  in  the  code  and,  when  we 
changed  the  sign,  good  agreement  was  found  with  Case  1  and  Case  6  as  shown  in 
Table  3.  However,  further  examination  revealed  the  sign  to  be  correct.  Although 
the  problem  could  not  be  pinpointed,  this  investigation  did  identify  the  model's 
weakness  in  the  IR  cooling  process. 

-L2  Planetary  Boundary  Layer  Depth  During  Stable  Conditions 

_i  22 

During  stable  conditions  (z  L  >  0),  Clarke  empirically  determined  the  depth 
of  the  PBL  to  be  approximately 

H  =  .  23u*/f  (7) 

where  u*  is  the  friction  velocity  in  cm  sec  1  and  f  is  the  coriolis  parameter  in 
sec  This  formulation  is  used  in  the  model.  However,  in  Case  4,  this  proved 
to  be  a  weak  point. 

In  Monin-Obukhov  theory  of  similarity,  a  scale  length,  L,  is  determined  based 
on  heat  and  momentum  fluxes  in  the  surface  boundary  layer.  In  this  model,  the 
surface  boundary  layer  extends  from  zq  to  2  m.  The  calculated  L  is  used  to  deter¬ 
mine  the  stability  of  the  layer.  Because  of  the  temperature  structure  in  the  lower 
2  m  for  Case  4,  the  model  calculated  a  stable  condition  and  used  Eq.  (7)  to  deter¬ 
mine  the  height  of  the  boundary  layer.  Based  on  u*,  the  boundary  layer  height  was 
determined  to  be  297  m  instead  of  500  to  560  m  as  specified  in  the  initial  and  verifi¬ 
cation  data.  The  free  atmosphere  is  assumed  to  exist  above  the  top  of  the  boundary 
layer  and  turbulent  intensity  is  assumed  to  approach  zero.  This  means  that  the 
cloud  layer  that  extended  from  420  m  to  500  m  underwent  no  turbulent  change.  As  a 
result  of  IR  radiational  cooling,  the  cloud  layer  warmed  slightly  and  evaporated  all 
the  liquid  water. 

It  is  interesting  to  note  here  that,  when  the  apparent  incorrect  minus  sign  in 
the  IR  radiational  cooling  routine  was  changed  (Section  4. 1),  the  results  for  Cases  1 
and  6  improved  greatly  but  those  for  Case  4  did  not.  The  cloud  was  not  evaporated, 
but  it  was  not  altered  (lowered  or  raised)  either.  It  continued  to  exist  between 
420  m  to  500  m  along  the  entire  horizontal  dimension  of  the  model.  This  was  caused 
by  the  boundary  layer  height  being  fixed  at  350  m.  Therefore,  the  cloud  existed  in 
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Table  1.  Description  and  Summary  of  Calspan's  Six  Data  Sets 


CASE 

TYPE 

FORECAST 

PERIOD 

(H) 

TIME 

OF 

DAY 

LOW  LEVEL 
TEMP 
STRUCTURE 

INITIAL 

CLOUD 

22  May  1978 
Case  i 

Stratus/ 
no  stratus 

4 

Morning 

Unstable 
to  Neutral 

yes 

2  Aug  1975 
Case  2 

Shallow 
cold  water 
advection  fog 

3 

Evening 

Inverted 

no 

7  Oct  1976 
Case  3 

Thinning 
stratus 
redeveloping 
to  form  fog 

10, 

14, 

20 

Morning 

to 

Night 

Neutral 

yes 

14-15  July 
1973 

Case  4 

Stratus 

thickening 

10 

Night 

Neutral 

yes 

5  Aug  1975 
Case  5 

Shallow  cold 
water  fog 
deepening  ove 
warm  water 

1-3 

r 

Morning 

Inverted  to 
Isothermal 

no 

29  Aug  1972 
Case  6 

Stratus 
lowering  to 
fog 

3 

Evening 

Neutral 

yes 

Table  2.  Comparison  of  Initial  and  Verification  Data  From  Calspan's  Six  Sets  With 
Model  Output 


r - .. - - -  -  -  ■  ~i 

« 

VERIFICATION /MODEL 

INITIAL 

HEIGHT 

TEMPERATURE 

DEPTH  OF  BOUNDARY 

j 

CASE 

CLOUD 

RANGE  (m) 

DIFFERENTIAL  (°C) 

LAVER  <m) 

In  it  lal 

Verification 

Model 

• 

* 

la 

no 

0-800 

8  to  -1.0 

790 

850 

780 

* 

lb 

yes 

0-800 

-0.  5  to  -1.  5 

790 

850 

863 

t 

2 

no 

0-28 

-3.  0  to  -6.  0 

Inverted 

inverted 

366 

i 

3.0 

yes 

0-250 

-2.0  to  3.0 

270 

250 

287 

3  14 

yes 

0-100 

-1.0  to  -2.0 

270 

100 

323 

320 

yes 

0-250 

-2.  0  to  -3.  0 

270 

250 

275 

4 

yes 

0-560 

0  to  +1. 0 

500 

560 

297 

51 

no 

0-28 

+3.  0  to  +4.  2 

inverted 

inverted 

177 

52 

no 

0-28 

0  to  +4.  0 

inverted 

17.  5 

190 

5- 

0-30 

0  to  +4. 0 

inverted 

30 

216 

*1 

3 

* 

cold 

yes 

0-200 

0  to  +2. 0 

200 

200 

161 

,1 

warm 

yes 

0-200 

0  to  4 2. 0 

200 

200 

161 

- 

» 

• 

26 
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Table  1  (contd).  Description  and  Summary  of  Calspan's  Six  Data  Sets 


r  -  ■  1  1  .  . .  -  - 1 

SEA  SURFACE  TEMP  | 

CASE 

WINDS 

SYNOPTIC  SCALE 
VERTICAL  MOTION 

Relative 
to  Air 

Trend  with 
Distance 

22  May  1978 
Case  i 

mod 

weak  upward 

warm 

warming 

2  Aug  1975 

Case  2 

mod 

cold 

cooling 

7  Oct  1976 

Case  3 

It 

strong 

subsidence 

warm 

warming 

14-15  July 

1973 

Case  4 

It 

zero 

warm 

constant 

5  Aug  1975 

Case  5 

mod 

cold 

to 

warm 

warming 

29  Aug  1972 

Case  6 

It 

zero 

warm/ 

cold 

Table  2  (contd).  Comparison  of  Initial  and  Verification  Data  From  Calspan's  Six 
Sets  With  Model  Output 


_ _ _ _ _ _ _  _ 

CASE 

FOG/CLOUD 

BASE  AND  TOP  (m) 

CLOUD/THICKNESS 

MAXIMUM  LIQUID  WATER  CONTENT 

Initial 
Base /Top 

Verification 

Base/Top 

Model 

Base /Top 

Thick 

(m) 

Initial 

Max  IAVC 
(g/m3) 

Verification 

Thick 

(m) 

Thick 

(m) 

Model 

Max  LWC 
(g/mty 

la 

688 

688 

468 

800 

None 

0 

0 

322 

0 

None 

lb 

464 

688 

468 

800 

None 

224 

.  36 

332 

0 

None 

2 

- 

- 

0 

26 

None 

0 

0 

28 

0 

None 

310 

144 

265 

None 

None 

121 

.20 

0 

0 

None 

314 

144 

265 

None 

None 

121 

.20 

0 

0 

None 

320 

144 

265 

44 

250 

None 

121 

.20 

208 

0 

None 

4 

420 

500 

250 

610 

None 

80 

.26 

300 

0 

None 

5i 

- 

- 

0 

50 

None 

0 

0 

50 

0 

None 
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the  free  atmosphere  and  underwent  little  change.  The  boundary  layer  height  deter¬ 
mination  overwhelmed  the  "improved"  IR  cooling  effects. 

4.3  Vertical  Motion 

In  the  model,  the  vertical  motion  is  derived  by  integrating  the  equation  of  con¬ 
tinuity  for  incompressible,  nondivergent  flow: 

w(z )  =  -J2  -^2.  dz.  (8) 

o  a  x 

This  is  a  definite  weakness  in  situations  such  as  Case  3,  where  strong  subsidence 
plays  an  important  role  in  the  evolution  of  the  atmosphere.  In  the  model  under  in¬ 
vestigation,  the  initial  subsidence  cannot  be  taken  into  account.  For  Case  3,  this 
has  a  significant  effect.  In  addition,  because  of  the  complete  dependence  on  du/dx, 
the  resultant  vertical  velocity  field  appears  almost  random. 

4.4  Short  Wave  Radiation 

In  Cases  1,  3,  and  5,  solar  radiation  effects  played  an  important  part  in  the 

development /dissipation  of  the  fog  and  clouds.  This  model  neglects  short  wave  ra- 

25 

diation  effects.  Stephens  et  al  reported  in  their  research  that  measured  short 
wave  flux  divergence  can  be  equal  in  magnitude  and  sometimes  larger  than  the 
long  wave  flux  divergence.  Therefore,  the  lack  of  solar  radiation  effects  is  a  sig¬ 
nificant  weakness. 

5.  SUGGESTIONS  AND  CONCLUSIONS 

The  results  from  the  numerical  model  appear  to  agree  very  well  with  the  re- 
15 

suits  reported  by  Barker.  This  agreement  indicates  that  the  model  realistically 
simulates  the  formation  of  fog  when  warm  moist  air  is  advected  over  cold  water. 
The  model  shows  that  a  stable  moist  layer,  when  cooled  from  below,  can  form  a 
surface  based  inversion  layer  that  remains  well  mixed. 

However,  as  described  in  Section  4,  a  minimum  of  four  changes  must  be  made 
before  the  model  can  be  used  outside  of  this  limited  context.  These  four  changes 
are: 

(a)  Include  solar  radiation  effects.  This  would  permit  diurnal  variations. 

(b)  Include  a  better  treatment  of  vertical  motions  and  a  way  to  input  initial 
vertical  motions.  This  would  permit  modeling  cases  where  subsidence  plays  a 
large  role. 

(c)  Develop  a  different  method  to  determine  the  height  of  boundary  layer  for 


stable  conditions.  This  task  would  be  easily  accomplished  and  would  produce  sig¬ 
nificant  benefits. 

(d)  Modify  the  infrared  radiative  transfer  equations.  This  is  the  most  impor¬ 
tant  change  of  the  four.  An  in-depth  review  and  necessary  changes  would  produce 
very  large  benefits. 

This  model  still  has  potential  in  an  Automated  Weather  Distribution  System 
(AWDS).  However,  the  weaknesses  described  in  Section  4  must  first  be  corrected. 
In  addition,  these  improvements  must  be  accomplished  keeping  in  mind  the  limited 
storage  and  processing  time  of  a  minicomputer. 
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Figure  A2.  Initial,  Verification,  and  Model  Temperature  and  Dew  Point  Temperature  Profiles  for 
Calspan  Case  1  at  4  Hours,  Cloud  Layer  Between  464  m  and  688  m 


Figure  A3.  Initial,  Verification,  and  Model  Temperature  and  Dew  Point  Temperature  Profiles  for 
Calspan  Case  2  at  3  Hours,  No  Initial  Cloud 


Figure  A 5.  Initial,  Verification,  and  Model  Temperature  and  Dew  Point  Temperature  Profiles  for 
Calspan  Case  3  at  14  Hours,  Initial  Cloud  Layer  Between  200  m  and  265  m 


Figure  A6.  Initial,  Verification,  and  Model  Temperature  and  Dew  Point  Temperature  Profiles  for 
Calspan  Case  3  at  20  Hours,  Initial  Cloud  Layer  Between  200  m  and  265  m 


Figure  A 9.  Initial,  Verification,  and  Model  Temperature  and  Dew  Point  Temperature  Profiles  for 
Calspan  Case  5  at  2  Hours,  No  Initial  Cloud 
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